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Abstract 

We consider, in the framework of General Relativity, the linear ap- 
proximation of the gravitational field of the Earth taking into account 
its mass, its quadrupole moment, its shape and its diurnal rotation. 

We conclude that in the frame of reference co-moving with the 
Earth the local anisotropy of the space is of the order of 1CT 12 - 1CT 13 
and could be observed. 



1 Introduction 

Three different types of experiments have been performed to test the isotropy 
of space. Although this concept has always a geometrical connotation what is 
really meant depends crucially on the type of experiments being considered. 

The first type of experiments are those of the type pioneered by Michelson 
and Morley ffl and developed by Miller 0, Joos || and many others using 
optical interferometers. Jaseja and al. || introduced a new setup using a 
stabilized laser as a standard of frequency and a Fabry-Perot as a standard 
of length, a technique that was used again by Brillet and Hall || who claim 
to have attained the best sensitivity. 

Historically the purpose of these experiments was to measure the absolute 
velocity of the Earth, i.e. the velocity with respect to that particular frame 
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of reference for which the light propagation, according to Newtonian Physics, 
was isotropic. In these experiments a(t) = Ac/c, where Ac is the difference 
of the round-trip speed of light along two perpendicular directions, say south 
and east, on a horizontal plane, is measured as a function of time t. These 
experiments [] have consistently yield results for max.a(t) ranging from 10~ 9 
to KT 13 . 

The second type of experiments was pioneered by Kennedy and Thorndike 
|| and is a variant of the Michelson-Morley one in the sense that light travels 
along a single direction fixed with respect to the Earth and one attempts to 
measure j3(t) = (c(t) — c(to)) / c(t ) , where c(t) is the round-trip speed of light 
along the fixed direction. 

The greatest value for a(t) obtained from the Michelson-Morley exper- 
iment is already 10 times smaller than the expected maximum Newtonian 
value 10~ 8 and therefore these experiments are considered to be good tests 
of Special Relativity. In fact, as Jaseja and al. have pointed out is not so 
obvious to decide what these experiments test. More precisely we ask here 
whether they should be analyzed in the framework of Special Relativity, ig- 
noring or not the centrifugal force due to the Earth's rotation?, or should 
this analysis have to take into account the Coriolis field and the Newtonian 
gravitational field in the framework of General Relativity? 

In the first case the Michelson-Morley experiment can be considered to 
be a test of isotropy of space in a cosmic global sense. In the second case 
the same experiment has to be considered as a test of the theory of light 
propagation in the presence of local inertial and gravitational fields. On the 
contrary the Kennedy-Thorndike experiment can only be considered as a test 
of the isotropy of space in the cosmic global sense. 

The distinction between the two points of view would be of paramount 
importance if a future experiment of one of these types yields a significant 
non-null result at some precision. The first point of view would mean the 
failure of Special Relativity while the second one does not mean that. It 
means that Special Relativity as a global model of space-time is not appli- 
cable here and that gravitational fields interfere with light already at a local 
level. 

The third type of experiments was pioneered by Drever and Hughes J7J. 

1 From a careful reading of |5j it follows that the upper limit is 10 -13 if one assumes 
that the anisotropy could have a local origin and 10~ 15 if one assumes that the anisotropy 
is necessarily a cosmic one. 
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They test the variation of an atomic or nuclear level splitting due to a fixed 
magnetic field as the orientation of this field with respect to the stars changes 
due to the Earth's rotation and motion around the Sun. These experiments 
yield a limit to this type of anisotropy of the order of 10~ 20 , a limit that 
Lamoreaux and al. f| pushed to 10~ 22 . Because the setup of these type 
of experiments is fixed with respect to the Earth they can be compared 
with the Kennedy-Thorndike type, although the physics which is involved 
is quite different. But they cannot be compared to the Michelson-Morley 
type because this one can discriminate between three cases: i) a(t) = 0, ii) 
a = constant 7^ and iii) a(t) a function of time. 

In this paper we analyze the local anisotropy of space in the neighborhood 
of the Earth's surface. The first section is based on the theory of Principal 
transformations. The analysis of the anisotropy of space is based on the 
exterior field only. The value of a predicted depends on the latitude of the 
site where the experiment is performed and three parameters close to the 
mean radius of the Earth. 

The analysis of the second section is global but model dependent. We have 
considered a model of the Earth with the observed oblateness but assumed 
its density to be constant. This analysis depends also on some heuristics that 
can be justified by the theory of Principal transformations. 

The third section contains a remainder of our interpretation of the aniso- 
tropy of space and a numerical comparison of the results derived from the 
approaches of the two preceding sections. 

2 The gravitational field of the Earth: A local 
approach. 

As an approximate model of the gravitational field of the Earth in the vacuum 
neighborhood of its surface we consider the following line-element 0: 



ds 2 = -(1 - 2U G )dt 2 + (1 + 2U G ){dr 2 + r 2 (dd 2 + sin 2 8d<p 2 )) (1) 



where U G is: 





2 We shall be using units such that c = 1 and G = 1. 
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M being the mass of the Earth, J 2 being the coefficient of its gravitational 
quadrupole moment and R the equatorial radius of the ellipsoidal surface. 
The "polar" coordinates (r, 9, <p) have been chosen for convenience and are 
related to a system of harmonic "Cartesian" ones (x % ) by the usual formulas: 

x 1 = r sin 9 cos (p, x 2 = r sin 9 sin (p, x 3 = rcos9 (3) 

Considering only those terms which are linear with respect to M and/or 
J2 the line-element [I] is a solution of Einstein's vacuum equations which 
actually reduce to the Newtonian equation: 

AU G = (4) 

To take into account the Earth's rotation Q we select a new frame of 
reference related to the preceding one by the transformation: 

ip -> ip + ttt (5) 

The line-element [1] becomes: 



ds 2 = -(1 -2(U G + U n ))dt 2 + 2A v dtd^ + (1 + 2U G ){dr 2 + r 2 (d9 2 + sin 2 6d<p 2 )) 

(6) 

where: 

U n = ^Q 2 r 2 sin 2 9, A v = Qr 2 sin 2 9 (7) 

Introducing new "cartesian coordinates" (x 1 ) connected with the new 
"polar" coordinates by the same formula |] the line-element ^ becomes: 

ds 2 = -(1 - 2U)dt 2 + 2A i dx i dt + (1 + 2U G )S ij dx i dx j (8) 

where: 

U = U G + U a , (9) 
and the non zero components of Ai are: 

A x = -Sly, A 2 = tlx. (10) 

From now on we shall neglect any terms proportional to powers of Q greater 
than 2. 
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Let r be the proper time along a time-like geodesic then, neglecting cubic 
and higher order terms in the velocities, the free fall of a test particle is 
described by the following equations: 

d 2 x k ~ k dx { dxi f dx'daA k k dxi 
-d^ + T ^^={ 1 + 6 ^^) A +2Q ^ 

where the Newtonian and Coriolis field are respectively: 

A* = ^d s {U G + U n ), = \5 ls (d s A 3 - d,A s ) (12) 

and: 

f * = SpiUa + StdjUa - 5, l3 5 ks d s U G - (Q^ + Q, k Q a )x l (13) 

Notice that, whether Q is taken into account or not, Eqs. |II| are meaningless 
without clarifying the meaning of the coordinates being used. Notice also 
that the newly defined "cartesian" coordinates x l are no longer harmonic 
coordinates. 

The symbols in |TT] are the Christoffel symbols corresponding to the 
Riemannian metric: 

ds 2 = ((1 + 2U G )Sij + n ik n jl x k x l )dx i dx j (14) 

which is the quotient of the space-time metric |6] by the Killing congruence of 
the frame of reference. 

This object is usually interpreted as defining the "geometry" of the space 
for the frame of reference being considered. This interpretation is correct or 
not depending on the meaning that one attaches to this word. It is certainly 
correct in a technical, current and accepted sense meaning a well established 
branch of mathematics. It is not correct if by "geometry of space" we mean 
that particular Riemannian geometry which allows the description of rigid 
objects which can be displaced, or compared when they occupy different po- 
sitions in space. This point has been strongly emphasized by Poincare with 
profound insight §|, and requires the geometry of space to have constant cur- 
vature [|I(J , the simplest case being of course the case in which this geometry 
is flat. 

If the object |n| is not the metric of space, what is it?. Before reminding 
a new interpretation which was first proposed in jy]] we are going to state 
and solve the problem of constructing the Principal transformation of the 
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Riemannian metric H, this being a concept introduced by one of us (Ll.B) 



in p|. 



Using polar coordinates we have: 

ds 2 = (1 + 2U G ) (dr 2 + r 2 d6 2 + r 2 sin 2 9(1 + 2U n )d<p 2 )) (15) 

The principal directions, i.e. the eigen- vectors of R\j, the Ricci tensor of 
this metric with respect to are to the required approximation proportional 
to: 

nii = (l,J 2 /,0), n 2i = (-Ja/,1,0), n M = (0,0,1) (16) 

where we have used the euclidean orthonormal co-basis (dr, rd8, r sin 9 dip) 
and where: 

R 2 ( r 5 \ 1 tt 2 R 3 

f = - 4 ^jy 1 -i^E) sindcosd with ^ = i^ ( 17 ) 

The vector rii is up to the required approximation, as a short calculation 
shows, collinear with the gradient of U pointing outwards. The vector n 2 
is therefore that vector lying in the meridian of the point being considered 
which is orthogonal to the vertical pointing south, and is the the vector 
lying in the corresponding parallel and pointing east. 

Using the preceding vectors as a basis of an orthogonal decomposition 
the line-element [15] can be written: 

ds 2 = (6 1 ) 2 + (6 2 ) 2 + (9 3 ) 2 (18) 

where: 



9 1 = (l+U G )(dr+J 2 frd6), 9 2 = (l+U G )(-J 2 fdr+rd9), 9 3 = (l+[/)rsin 

(19) 

By definition, in this particular Principal transform of [18] is a metric 

which can be written as: 



ds 2 = c 2 1 (9 1 ) 2 + c 2 2 (9 2 ) 2 + c 2 3 (9 3 ) 2 

the three functions Ci(x^) being chosen such that: 
i) the metric EH is euclidean, i.e. such that: 



(20) 
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ii) the following "quo-harmonicity condition" is satisfied: 



(21) 



(22) 



where T) k are the Christoffel symbols of the line-element |20], and g^ k is the 



inverse matrix of gjk- 

Notice that both groups of conditions are tensor equations and therefore 
the solutions Cj(x J ) are scalar functions independent of the system of coordi- 
nates being used, as long as this system is adapted to the Killing congruence 
defining the static frame of reference. 

Integrating pi] and 1221 we have obtained to the required approximation: 
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where i?i and R 2 are in this approach two model-dependent free parameters. 
They should be determined by matching the preceding results with the cor- 
responding ones in the interior of the Earth. This would require a global 
C 2 model including the interior gravitational field and therefore also a C° 
behavior of the density of the Earth. To deal with a discontinuous density 
requires to adopt a different approach. This is done in the next section. 
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3 The gravitational field of the Earth: A global 
approach. 

To begin with we consider in this section a model of the Earth with constant 
density taking into account its oblateness. More precisely we shall assume 
that the surface of the Earth £ is r 2 = a 2 (l + A 2 sin 2 6), where a is the polar 
radius, and A = \/ R 2 — a 2 / R the eccentricity, R as before is the equatorial 
radius. 

An elementary integration of the corresponding Poisson equation yields 
up to order J 2 the Newtonian potential: 



Ur 



— 1 + — — (1 - 3 cos 2 9) ) outside 
r \ 2 r 2 



i£ + Sliced ]] inside 

R\2 2R 2 2 1 2 2 R 2 R 2 



(23) 



which is of class C 1 across S, where J 2 = A 2 /5. 

Taking into account also the diurnal rotation of the Earth, the gravita- 
tional field in a comoving frame of reference is then described by the line 
element |], where now Uq is given by |23] instead of § and correspondingly we 



have to make the same substitution in [15. 

On the surface of the Earth E the Ricci tensor of the metric [15] is dis- 
continuous as it is the density and therefore the principal directions derived 
from the exterior field do not coincide with those derived from the interior 
one. This prevents the existence for this metric of a global principal trans- 
formation continuous across S. 

When a principal transformation exists then from E2I it follows that carte- 
sian coordinates of [2^, i. e. those with Y l - k = are quo-harmonic coordinates 



of [15], i. e. they are such that T l j k g lJ = 0. The approach will then be 
heuristically justified by this result. We shall obtain a global system of quo- 
harmonic coordinates, unique up to euclidean transformations, and we shall 
consider the anisotropy of the space to be the anisotropy of [l^ with respect 
to the euclidean metric which in this system of coordinates has components 
9ij = Sjj^ 

3 This point of view was already used in JlJ] where harmonic coordinates of the space- 
time were considered instead of the quo-harmonic coordinates of the space. It could be 
considered also in connection with rotating Fermi coordinates based on the world-line of 
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Under an infinitesimal coordinate transformation 

x* = z i + C(z j ) 

where Q are of the order appropriate to the approximation that we consider, 
the Riemannian metric |14] becomes: 



^, (lV )) + § + f§ 



The new coordinates Zi will be quo-harmonic if: JL3 



A ^) = f^ ( 24 ) 

We can solve this equation splitting the functions (k in their gravitational 
and rotational part (k = (j? + Ck- Then 

ACf=ff, (25) 

= -2Q kj Q j lZ l . (26) 

The solution for (f? can found by integrating the Poisson equation ^ 
because the "density" dUcldz k is known. Requiring the gradient of to be 
zero at infinity. In the new polar coordinates, i.e. coordinates (r, 9, ip) such 
that: 

z 1 = r sin # cosy?, z 2 = r sin 9 sin ip, z 3 = r cos 9 (27) 



we obtain: 
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the center of the Earth. These two coordinate dependent models are heuristic as the one 
presented in this section. The stronger legitimacy of the latter comes from being ackin 
with the method of the first section. 
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The rotational part can be obtained very easily from an inspection of 
the explicit expression of |56. Demanding the appropriate symmetry and the 
correct limit on the axis of symmetry we obtain: 
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(30) 



Writing the metric 13 as g+j = 5ij + hy and using polar coordinates we obtain 
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(31) 

Making more precise what we said before we define the eigen-values and 
the principal directions of the anisotropy of the metric [15] as the scalars c a 
and the vectors n % a solutions of the algebraic equations: 



(9ij ~ cJi^nl = 
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with hij given above by |3T]. To the required approximation the scalars c a are: 

c x = l + h rr) c 2 = l + ^, c 3 = l+ 2 V 2fl (32) 

r z sin # 

and the components of the principal directions in the euclidean orthonormal 
co-basis (dr, rd8, r sin 8d(p) are proportional to: 

ni i = (l,/(r,0),O), n 2i = (-/M),l,0), 71* = (0,0,1) (33) 

where: 

/M) 



M(l -3i? 2 /(5r 2 ))' 

Conclusion 



In section 2 we exhibited the anisotropy of the riemannian metric [15] with 
respect to an intrinsically associated euclidean metric. This association de- 
pends on three model dependent parameters Ri, R 2 and R. If the quadrupole 
moment were zero the three parameters would coincide and assuming a model 
of the Earth with constant density their value would be its radius. Therefore 
we can safely guess that the three parameters will be close to each other and 
of the order of the mean radius of the Earth. It is worthwhile to notice that 
on any horizontal plane, i.e orthogonal to the gradient of U the anisotropy 
depends only on R. And in particular if r = R « Ri R 2 we get: 

a = c 2 - c 3 = ~ " 14^ 2 i? 2 ) sin 2 9 (34) 

Assuming the following numerical values: M = 0.00444 R = 6378164 and 
Q = 2.434 x 10 -13 , the three of them measured in meters, and J 2 = 0.0010826 
dimensionless, at a latitude of 45° we obtain: 

a = C2 - c 3 = 2.5 x 10~ 12 (35) 

In section 3 we constructed a globally C 1 model of the gravitational field 
of the Earth assuming it to be of constant density and having an ellipsoidal 
shape.0 Then we obtained a global system of quo-harmonic coordinates of 

4 This model of the Earth would yield a value of A 2 = 0.00541. while the observed value 
is A 2 = 0.00674 
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the space. It follows from that that we can associate with the global field 
of the Earth an euclidean metric with components = 5^ in this unique 
system of quo-harmonic coordinates. This is again a model dependent but 
intrinsic association. 

With this model the anisotropy in the horizontal plane, for r = R is: 

4MJ 2 , 3 o2d2 \ . , 



a = c 2 - c 3 = [-j-jf + 2 Q R ) sin 6 ( 36 ) 
and with J 2 , M, R, Q and 6 as before, this formula yields: 

a = c 2 - c 3 = 6.2 x KT 13 (37) 
The two results in |35| and [37] differ, although by less of a factor of 10. 



The safe conclusion that we can draw from our precedent analysis is that it 
predicts an anisotropy on the horizontal plane in the range 10~ 12 — 10~ 13 . f] It 
is however important to remind that the theory of Principal transformations 
tell us that for a model of the gravitational field of the Earth of class C 2 the 
two methods would yield the same result. 

The important question is now the following: how could this anisotropy 
show up?. Because the anisotropy is not due to any cosmic effect but is 
due to the gravitational field of the Earth and to its diurnal rotation and 
because it is relative to a co-moving frame of reference, no experiment of 
the Kennedy-Thorndike or Hughes-Drever type could reveal it. If instead, 
as it was first proposed in ||11|| , and [113], we renounce to interpret the 



metric [T| with components g^ as the geometry of space, a role which belongs 
to the euclidean metric g~ij, and we interpret the eigen- values of the first of 
this objects with respect to the second as the velocities of light along the 
principal directions of an anisotropic propagation, then the experiments of 
the Michelson-Morley type should reveal such anisotropy. 

As a matter of fact Brillet and Hall, the authors of the most recent repe- 
tition of this experiment obtained a raw result of the order of 2.1 x 10~ 13 at 
the latitude of +40°. If confirmed our model could provide an explanation 
for this result. In our opinion is an urgent task that somebody repeats this 
experiment using their setup or any other more appropriate variant. 



The harmonic coordinate method and the Fermi coordinate one yield values in the 
same range. 
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